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We study dissipative effects due to inertial forces acting on matter fields confined to accelerated 
boundaries in 1 + 1, 2 + 1, and 3 + 1 dimensions. These matter fields describe the internal degrees 
of freedom of 'mirrors' and impose, on the surfaces where they are defined, boundary conditions on 
a fluctuating 'vacuum' field. We construct different models, involving either scalar or Dirac matter 
fields coupled to a vacuum scalar field, and use effective action techniques to calculate the strength 
of dissipation. In the case of massless Dirac fields, the results could be used to describe the inertial 
forces on an accelerated graphene sheet. 
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I. INTRODUCTION 



In what is one of the most startling manifestations of 
the quantum nature of the electromagnetic (EM) field, 
the Dynamical Casimir Effect (DCE) predicts that, in 
the presence of an accelerated mirror, the quantum vac- 
uum evolves to an excited state (i.e., one with a non- 
vanishing number of photons). As a consequence, the 
mirror will experience a dissipative force. Both the dis- 
sipation and the number of created photons are rather 
small for a single accelerated mirror. However, the effect 
can be enhanced in electromagnetic cavities with time 
dependent length or, more generally, with time depen- 
dent electromagnetic properties (see the reviews [l|, and 
references therein). 

In this letter we consider an effect that, although of 
a different nature than the DCE, may nevertheless arise 
under fairly similar conditions: acceleration of the mir- 
rors. The main difference with the DCE is that we shall 
consider the 'inertial' dissipation that results from the 
transfer of energy to the matter fields inside the mirror, 
rather than to the EM field. 

Effects due to the inertial forces in accelerated me- 
dia are known, and usually described in terms of in- 
duced time-dependent polarizations or currents, like in 
the Barnett and Stewart-Tolman effects [2| ■ More gener- 
ally, there will be inertial forces as long as there is energy 
transfer from the collective coordinates of the mirrors, 
that specify their position and shape as a function of 
time, to their internal degrees of freedom. 

Contrary to the vacuum friction effects considered up 
to now [H-[l]) the inertial effects described in this letter 
will exist even when the internal degrees of freedom are 
not coupled to the EM field. Our aim in this work is to 
evaluate these effects using specific models for the matter 
fields, allowing for mirrors with a quite general space and 
time dependence. 



II. THE MODEL 

Let us begin by considering a scalar model, described 
by an Euclidean action S with the structure: 

S = S f + S m + Si , (1) 

where Sf describes the vacuum field dynamics in the ab- 
sence of mirrors. It will be assumed to be a free field 
theory in d+ 1 dimensions (d = 1, 2 or 3). S m is the free 
part of the matter field action, i.e., it ignores the coupling 
to the vacuum field; this coupling is instead contained in 
the last term, Si. 

Both S m and Si have support on AA., the d-dimensional 
spacetime manifold swept by S(r), the space occupied 
by the mirrors at the time t, in the course of their time 
evolution. 

The vacuum field if is assumed to be a scalar, equipped 
with the action 

SfQp) = \J d d+1 x dM*)dMx) • (2) 

In this scalar model, the matter degrees of freedom are 
described by £, also a real scalar field, but living on 
AA, a curved d-dimensional manifold. To make this 
more explicit, we introduce d variables: cr°, a 1 , ... <r d ^ 1 , 
such that a is the time coordinate; AA can then be 
parametrized by d + 1 functions /i = 0, 1, . . . , d, as 
follows: 

M: a — > j^V,....^- 1 ) (3) 

such that y° = <r° = x° = r and y- 7 = y^(r, a 1 , . . . , a x ) 
(j = 1, . . . d), is a parametrization of S(t) (for each r). 
Note that AA may denote either one or more than one 
mirror. When dealing with more than one mirror, we 
shall introduce an alternative, more explicit notation, 
obtained by using a different parametrization for each 
connected piece in AA. 

We assume the scalar matter field action to be given 
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by: 

SmfcM) = \J d d ay/gja)[g^{a)d a ^a)d^{a) 

+ few] , (4) 

where g a p{&) is the induced metric on A4: 

g a0 (a) = d a y»(a)d f} y»(a) , a = 0, 1, . . . , d - 1 , (5) 

g{(j) its determinant, and fi the mass of the field. 

For the Si piece, on the other hand, we shall consider 
a term of the form: 

Sj(^tM) = -<C / <?<rViiFiZ{<r)<p[v(p)] > (6) 



where £ denotes the coupling constant. 

The action 5 can be considered as a toy model for the 
interaction between the EM field (ip in the toy model) 
and charged degrees of freedom living on the mirrors (de- 
scribed by £), although, as it shall be evident below, it 
is unnecessary to have a coupling between the internal 
degrees of freedom and the vacuum field for dissipation 
to arise. It is nevertheless interesting to consider models 
where that coupling exists, in order to compare 'inertial' 
and DCE dissipations. 

The scalar nature of the matter field is appropriate to 
the field configuration describing an 'internal' (charge) 
degree of freedom. We shall also consider the case of a 
Dirac-like matter field, but we shall not dwell with vector- 
like matter fields, which could account for phonon- like- 
excitations, and certainly also exhibit inertial dissipation. 

To take into account the effects of the internal degrees 
of freedom, we introduce T m (cp), the effective action due 
to the matter fields, defined as follows: 

e -r m ( v -M) = J £>£ e -S m (t,M)-Si(v,Z;M) ) (7) 

which also encompasses the dependence on the geometry 
oiM. 

Because of their quite different properties, we split T m 
into two contributions, as follows: 



where 



r m (cp;M) = Ti(M) + r b (<p;M) (8) 



Ti(M) =T m (ip;M)\ v=0 , 



T b (<p;M) =T m (ip;M) - T m (<p; M)\ v=0 . 

Each one of these terms is relevant to different effects: 
Ti(M) can induce dissipative effects, while T b (<p;M) is 
responsible for the emergence of boundary conditions on 
ip. Our focus shall be the calculation of Ti(A4), although 
we shall also derive the form of T b (ip; M) for the sake of 
completeness. 



It is quite straightforward to see that: 

I 



Ti{M) 
where JC is the operator 



• Tr In K 



(9) 



K 



d a \g 1/2 g ap dp 



,1/2 



(10) 

and Aju = g~ 1 / 2 d a g 1 ^ 2 g a ^ l dp is the Laplacian cor- 
responding to the induced metric on M.. Note that the 
operator K. is Hermitian for the scalar product (/i, /a) = 
jd d a(h{a))*h{a). 

Besides, the boundary interaction term Y b {ip;M), is 
given by: 

T b (^M) = £ J d d ag 1/2 (v) J d d a'g^\a') 

x ip[y{aj\lC- l {*,a')<p[ y {a'j\ . (11) 

Expressions Q and (|lip . properly interpreted, are 
valid for one or more than one mirrors. If there is just 
one mirror, their interpretation is immediate: one needs 
to know a parametrization of its surface, then all the ob- 
jects involved in © and (|TTj) are derived from it and the 
free scalar field propagator. When there are N > 1 mir- 
rors, we introduce a discrete index a, and use a different 
notation for the functions that parametrize each mirror: 
y%, a = 1, . . . , N (which can be, for example, restrictions 
of the parametrization of Ai) . Then, since the action for 
the matter fields is localized on each mirror, we see that 
the effective action becomes the sum of the contributions 
due to each mirror: 



N 



N 



Ti(M) = r l (M a ) , r b ( w m) = r &(^; M a ) , 

a— 1 a— 1 

where M. a is the world- volume of each mirror. 

In other words, quantum effects of the matter fields 
are, for this model, additive with respect to the mirrors. 
Thus, in the rest of this paper we shall only calculate I\ 
and T b for single mirrors. Before presenting results of the 
evaluation of those contributions for different numbers of 
dimensions, we check that, in the limit when the micro- 
scopic degrees of freedom do not propagate, the bound- 
ary interaction term reduces to the standard one. Indeed, 
keeping just the leading term in a large- expansion: 



M - 2 g-^{a)S{a-a') 



(13) 



so that for the leading behaviour of the boundary inter- 
action term we have the following expression: 



X 



dW\v)v 2 Wj\ (14) 



where A = . This is the usual 5-like potential that 
has been used to analyze the Casimir effect with imper- 



fect boundary conditions [5|. 
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III. NO INERTIAL DISSIPATION IN 1 + 1 
DIMENSIONS 

To gain some insight into the previous expressions, let 
us first evaluate the effective action for a case where the 
outcome will be that there will be no inertial dissipation, 
namely, d — 1 . Here the vacuum field is a massless real 
scalar field (in f + f dimensions) while £ is a single quan- 
tum mechanical degree of freedom. For the model we 
consider, £ cannot experience inertial forces, since it is 
not spatial (its values belong to an internal space), and 
moreover it is only a function of time. 

The spacetime manifold M. is one-dimensional, a plane 
curve which we parametrize as follows: r — > (r,q(r)), 
where r is the time and q(r) the position of the mirror. 
The S m and Si terms may therefore be written in the 
form: 



•S,,ii':.V!: = ~ / ,Jr 



V)(^) 2 + M 2 e( T)eV) 



Si(<p,S;M) 



i( / d,Te{T)£,(T)(p(T,q(T)) 



(15) 



where e(r) = \/g(r), and g[r) = 1 + g 2 (t) . Then we 
note that, by performing the reparametrization: r s, 
such that ds = e(r)<ir, we get: 



-Ids 



(^) 2 +. 2 ?(s) 



i( / ds£(s)ip(s,q(s)) 



Sj(<p,S;M) 

where we have introduced the notations: 

£00 = €[t(s)] , <p[s, q(s)] = (p[r(s),q(s)] , 



(16) 



(17) 



Integrating out the matter field (we assume its path in- 
tegral to be invariant under reparametrizations |12j). we 
see that in the effective action term 1^, the dependence 
on e(r) (hence on q(r)) is completely erased. Thus no 
dissipative effects from this origin may arise in d = 1, as 
advanced. 

Regarding the Tf, term, we find that: 



e 

2 



ds 



d S 'v(s,q(s))A( S - S ')v( S ',q( S ')) 

(18) 

with A(s — s') = e ^ s s '/2/i. This boundary interac- 
tion term can be also written in a Fourier representation, 
so that it produces a frequency dependent coupling be- 
tween the vacuum field and the mirror: 



e 

2 



du> 



1 



2tt uj 2 



(19) 



where J- denotes Fourier transform with respect to the 
s variable. We note that the outcome of this study 
amounts to the property that the Tjj term is, essentially, 
a reparametrization invariant interaction, which is tan- 
tamount to using a relativistic invariant term. For the 



non-dynamical limit, this is precisely one of the cases we 
have studied in [6|. Indeed, when the matter degrees of 
freedom do not propagate, 

n(y,M) ~ \ J d s [^( s ,q(s))] 2 
A 



drVl + g 2 (T)^(r,< ? (r))] z ,(20) 



with A = which is identical to equation (33) of [f|. 

In the general, nonlocal case, we can still find the ef- 
fective action r[g(r)], for a single mirror, obtained by 
performing now the functional integration over the vac- 
uum field: 



e -r[?(r)] = / VLpe -S s ( V )~T b { V ;M) 



(21) 



To second order in the mirror departure from its average 
position, we get: 



dre(r) / dT'e(r')q(r)F s (r-T')q(T') 

-oo J — oo 

(22) 



where: 

F s (t-t') 

and: 



du> 
2^ 



e i*(T-T>) Fs ^) , Fs ^) = F(w)-F(0), 

(23) 



1 



V ' 4 J 2ir 1 C 2|;/ + w| J 1 1 

(24) 

Although the integral cannot be evaluated exactly, it is 
clear that its UV behaviour is much improved with re- 
spect to the one corresponding to a local Tf,, like in [f|. 
The physical reason for that is that the nonlocal Tj, intro- 
duces a kind of cutoff frequency (of order fx) above which 
the mirror becomes ineffective to impose the boundary 
conditions. 



IV. 2 + 1 DIMENSIONS: INERTIAL 
DISSIPATION ON A MOVING STRING 

The situation is more interesting regarding dissipation 
when d = 2, since the space dependence of the matter 
field allow for the action of inertial forces. Besides, if 
fi = 0, it is possible to derive a closed-form expression 
for T m (Ai). Indeed, this object is formally identical to 
the effective action for a massless real scalar field (£) 
in two dimensions, which can be found exactly, and it is 
nontrivial only if R, the scalar curvature of M. is different 
from 0. Thus, before writing the explicit form of the 
result, we find conditions for R to be different from 0. 

Using y to denote the two spatial components of the 
parametric form of .M, namely, y = y(<T°,cr 1 ) (cr° = r, 
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a 1 = a), we find, after a quite straightforward calcula- 
tion, that: 



R 



A 



(25) 



where g is the determinant of the induced metric, and: 

A = (d iy rd 2 y 2 d 2 iy 2 ~ dijMV c&My 2 

- d'y 1 8 1 y 1 8 iy 2 Ofy 2 + (d lV 2 ) 2 3 V d^y 1 

- (d m ) 2 (d d iy 2 ) 2 - (d iy 2 ) 2 (dod^ 1 ) 2 

+ 2d 1 y 1 d lV 2 dod.y 1 d 8 iy 2 , (26) 

where y 1 and y 2 are the two spatial components of y. A 
simpler expression for R may be written when the surface 
can be written as a Monge patch (which is then neces- 
sarily open), that is, when A4 can be parametrized as: 



y 



t , y 



a , y 



V(r, a) 



(27) 



Note that y(r, a) describes a moving open string in 2 + 1 
dimensions. In this case we have 



R 



F 



(28) 



with 

F = 



d 2 y d T d a y 
d a d T y diy 



, g=l + {d T y) 2 + {d a y? . (29) 



A necessary and sufficient condition for F to vanish can 
be obtained: 

F = <-> 3a E E / (d T + ad a )y(r,a) = . (30) 

This includes (since a can also be or oo), the cases of 
an F which depends only on r or only on a. We may 
write then all solutions to the condition above, for R to 
be zero (and T m trivial) as follows: 



V = f(r-a V) , 



(31) 



(for any a). Thus, there will be no effect if the evolution 
of the mirror has the form of a packet that evolves undis- 
turbed, with any speed (eventually zero). That is to say, 
they are solutions of the wave equation with wave speed 
a (arbitrary), and with a definite chirality: namely, they 
should either move to the left or to the right, but cannot 
be a combination of both. The absence of dissipative ef- 
fects for this kind of evolutions is certainly related to the 
Lorentz invariance of the model. 

On the other hand, combining both left and right 
movers, one certainly can have an F =/= 0. Thus, what 
is perhaps the simplest example of a configuration where 
F =/= corresponds to standing waves. 

We now recall the form of the (exact) effective action 
for a massless real scalar field in 1 + 1 dimensions [7j , to 
write: 

T m (M) = f d 2 xR^R. (32) 



We shall, in what follows, only consider the massless mat- 
ter field case. The reason is that for a massive field the 
corresponding contribution is strongly suppressed, unless 
the curvature is big in comparison with a scale set by the 
mass /i. In such a case, however, the massless field case 
should be a good approximation. Moreover, in the mas- 
sive case there are also local terms in the curvature, but 
their locality makes them irrelevant to dissipative effects. 

Of course, inserting the expression corresponding to a 
Monge Patch, one gets a more explicit form. Doing that, 
one may look for the T m that results from the lowest 
order contribution, in a derivative expansion of y(r,a): 

T TO ~ J drda J dr' da' F (t , a) I^ 1 (t - t' ,<t - a') 
x F(r',cr'). (33) 

where A^" 1 is the free propagator (which has the usual 
logarithmic form), and F has been defined in 

The pole in the free propagator is responsible for the 
existence of a nontrivial dissipation. Indeed, we may 
Fourier transform, and rotate to real time, to obtain the 
imaginary part of r m : 

ImTrj ~ — / — r^-r \\F(k = |fci|,ki)| 2 
L J 96tt J 2ir fcJ L 1 v 1 11 



\F(k 



(34) 



where F(ko, fci) is the Fourier transform of F. Note that 
the imaginary part sees the arguments of this function 
on-shell. On the other hand, F is quadratic in y and 
therefore the dissipation is quartic in the deformation of 
the string. 



3 + 1 DIMENSIONS: INERTIAL DISSIPATION 
ON A MOVING MIRROR 



We conclude our study of the scalar field model with 
the d = 3 case. Here, it is not possible to find an exact 
effective action for an arbitrary manifold M. (even when 
fi = 0). The effective action can, however, be computed 
approximately, in an expansion in powers of the curva- 
ture. One has, up to the second order Q: 



1 



64 x 2 3 / 2 
+ a 2 R(~A)-iR +r ioc „/ . 



d 3 ay/g(a) ai Rap (- A) * R a p 

(35) 



where a\ = — 1 and a 2 — 1/8. Ti oca i is a divergent con- 
tribution constructed with the induced metric and its 
derivatives, and we will assume that these divergences 
are absorbed into appropriate counterterms in the clas- 
sical action. Since local terms are not relevant for the 
dissipative effects considered here, we will omit Ti oca i in 
what follows. 
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The nonlocal kernel (—A) 2 can be formally defined 
using the integral representation 



dm 



1 



TO 



(36) 



in terms of the two-point function of a massive scalar 
field (-A + m 2 )- 1 . 

Consistently with the expansion of the effective action 
in powers of the curvature, we shall analyze the nonlo- 
cal part of Ti(A4) under the assumption that the mir- 
ror is almost flat, and that the time dependence of its 
surface is smooth. This means small departures from a 
flat surface almost flat manifold M.. To implement this 
in a quite straightforward way, we first assume that the 
parametrization used is of the form: 



2/(<xW 2 ) 



for jji 
for /j, - 



0,1,2 
3 . 



(37) 



Then it follows that g a p = 5 Q( 3+/i Q /3, with h a p = d a ydpy. 
We then perform an expansion in the deformation of the 
manifold, i.e., around the 5 a p metric: g a p = 5 a p + /i Q( g, 
assuming that h a p is small. To find the lowest non-trivial 
contribution in an expansion in h, we note that: 



^372 / d 3 <rR aP (-A yi 



64 x 



(38) 



where R a p is the Ricci tensor expanded to its lowest non- 
trivial order in h a p and (— A )~^ is the nonlocal kernel 
for a flat spacetime. 

To linear order in h a a we have 



Using this equation, it is straightforward to write the 
Ricci tensor in terms of the function y(a). For instance 
we have 

Rqo ~ d diy d diy + d d 2 y d d 2 y - d d y d^iy 
- d d yd 2 d 2 y , (40) 

and similar expressions for the other components. 

As in the two-dimensional case, one can readily show 
that all components of the Ricci tensor vanish when 
the evolution of the mirror is described by a packet 
that moves with constant velocity in any direction, i.e. 
y = f(a° — oT x a % ). Therefore, the manifold M is flat 
for this kind of evolutions, since in d = 3 the Riemann 
tensor is entirely determined by the Ricci tensor. Once 
more, simple evolutions that produce dissipative effects 
are standing waves. 

The dissipative effects induced by the matter degrees of 
freedom are related to the imaginary part of the effective 
action in Minkowski spacetime. Rotating to real time, 



and introducing Fourier transforms for the fiat propaga- 
tor and the Ricci tensor we obtain 



Im [Fi] 



Illl 



1 



dm 



d 3 k 



32tt2 3 / 2 J q ""'"J (2tt 
Ra^Rpai-k^rfPr]^ 



( A: 2 , + k\ + k 2 



to 



(.41) 



Introducing the notation 



\B(k)\ 2 = R aP {k)R pa {-k)[ifP^ - -t^tH , (42) 
and performing the integration over ko, we find 



Im [r< 



1 



dm 



d 2 k 



2 3 / 2 64 7 J (2tt) 2 
|B(fc = Wk ,k)| 2 + |B(fc = - Wk! k^ 2 



{43) 



with «k = \/k\ + k\ + to 2 . The structure of this result 
is similar to that of the two-dimensional case (1341) , except 
for the additional integral in the parameter to, used to 
have a suitable integral representation of the nonlocal 
kernel (— A 2 )~2. 

It is worth mentioning that a case in which the matter 
field is a massless Dirac field ip (in 2 + 1 dimensions) leads 
to very similar results. Assuming the matrices 7 (a = 
1, 2, 3) for this field to be in a reducible representation of 
Clifford's algebra, chosen in such a way that the parity 
anomaly is cancelled, the matter field action S m is given 
by: 



S m (^,^;M) = d 3 ay/g{aj^a)D^(a) (44) 



where D = 7 a e^(cr)D At , is the dreibein and = 
(39) d, I^b^ab- Here, lu^ denotes the spin connection 



(determined by the dreibein and Christoffel symbol) and 
(T a t, = 5 [7a j 7b] • The result for Ti may also be read 
from [8j, and it is given by (|35j) with the only differ- 
ence being in the numerical values of the coefficients a\ 
and a 2 , that we will not need in what follows. Thus, all 
the estimates below, as well as the conditions for no dis- 
sipation hold true for this model as well. Note that the 
Dirac field model mig ht be relevant to models describing 
graphene surfaces [lOf . 

We can obtain an estimation of Im [r^] for the scalar 
and Dirac fields, by assuming that the dynamics of the 
mirror is described by a standing wave of the form 



y(crVV 2 ) = ?/o cosher ) cos (a 1 j L) 



(45) 



where we have denoted by yo the amplitude of the oscilla- 
tions, O its frequency and L the typical distance between 
nodes in the direction of a 1 . For simplicity, consider evo- 
lutions of the mirror that are invariant along a 2 . From 
and ((10]), we can calculate the Ricci tensor. For 
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example, we can estimate the Fourier transform of (1401) . 
which can be written as 



4 3 2q2 

R 00 (k) = f^— S(k 2 )[S(h)(S(k 

+ S(ko + 2n)) + 5(k Q )(6(k 1 -2/L)^ 



- 2ft) 

5{kx+2/L))]. 



Similar expressions can be found for the 11-component of 
the Ricci tensor, and it is easy to see that R\o = i?oi = 0. 

In order to have a rough estimation of the magnitude 
of Im [l\] when the mirror's dynamics is approximated 
by Eq. (|45p . we need to perform integrations in Eg. (1431) . 
After integrating in k± and k^, we get 



Im [Ft] 



■yd 5 



L 4 



(46) 



where T is the total time during which the mirror is mov- 
ing, and £ denotes the total surface of the (unperturbed) 
mirror. The presence of these coefficients is a byproduct 
of having considered a delocalized standing wave as the 
one in f|45[) . It is worth noting that, for this particular 
motion of the mirror, there is no threshold for the dissi- 
pation effects. 

Restoring h and the propagation velocity (vf) of the 
matter fields inside the layer, the estimated magnitude 
of Im [rA becomes 



Im [Ti] hy^Vl 3 



(47) 



We may compare this result with the one obtained for 
one perfect mirror oscillating according with Eq. (|45|) . 
In this case, the DCE produces an imaginary part in the 
effective action as long as fl > c/L fllj ]: 



This ratio may be bigger or smaller than one depending 
on the value of the different parameters. It is interesting 
to remark that while the inertial dissipation is quartic in 
the amplitude yo, the dissipation associated to the DCE 
is quadratic. Moreover, for the particular motion (J45j , 
the inertial effects do not have a frequency threshold, 
while the DCE does. 

In Eq. (|49[) we compared the inertial dissipation with 
the DCE for a perfect mirror. Note, however, that, being 
independent of the coupling to the vacuum field, the in- 
ertial effect may be more relevant than the DCE in cases 
where that coupling is rather weak, namely, when the 
mirrors are far from perfect. 

It is also worth noting that the result mentioned above, 
that a deformation traveling undistorted (with a constant 
direction) does not produce inertial dissipation, may be 
relevant to the DCE. Indeed, in the case of having more 
than one mirror with such deformation, the configuration 
is free of inertial dissipation, while the DCE does not 
vanish. 

Finally, we wish to point out that we have considered 
here the effect due to inertial forces just for one matter 
field; in a real medium many different contributions can 
appear, coming from the excitation of different internal 
degrees of freedom. All of them will contribute to the 
imaginary part, with different parameters and kinemati- 
cal factors. 



im[r DCE ] ^ (1 __ § ! f)5 /2 



(48) 



When fl is above threshold, the ratio between both effects 
is approximately given by 
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